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ABSTRACT 


The performance characteristics of a residential solar 
assisted water-to-air heat pump heating system are analyzed 
for State College, Pennsylvania. A realistic residence and 
solar assisted water-to-air heat pump system are modeled for 
this northern climate using’ the transient simulation 
computer code TRNSYS developed by the University of 
Wisconsin- The system is studied over a one month winter 
period, December, using actual hourly weather data for State 
College. The system is analyzed for both the cloudiest and 
clearest December weather recorded in the last 3Q@ years. 
The collector area and storage tank capacity are varied and 
the effects on system performance are _ studied. Other 
parameters varied and studied are the collector/ storage 


tank heat exchanger effectiveness and the room temperature. 





CHAPTER 1 


a Introduction 


With the energy crisis worsening and conventional 
heating fuels becoming scarce and expensive a need exists 
for the development. of alternative heating fuels and 
systems. Solar energy which is "free" and readily available 
is a promising source of building heat. Unfortunately for 
northern climates direct solar heating of buildings can 
provide’ only a small portion of the total heating 
requirements without extremely large collector areas’ and 
heat storage capacity[1l])-.* The solar energy can however be 
used to augment conventional heating systems or used _ to 
improve the performance of the system, as with a heat pump. 
The heating capacity of the heat pump decreases as- the 
evaporator temperature (1-.e- the outdoor ambient air 
temperature) decreases. At the same time the heating load 
of the residence increases. The point at which the heat 
pump capacity equals the heating load is called the balance 
point. The heat pump is an economical heat source as long 
as it can operate at a temperature above its balance point. 
In cold northern climates the evaporator inlet temperature 
of the standard air-to-air heat pump installation is often 


below the balance point. For water-to-air heat pumps an 


* Numbers in brackets refer to references listed at the end 


of this paper. 
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existing warm water source (temperature above’ the balance 
point) for the heat pump evaporator is difficult to find in 
northern climates particularly for residential 
installations. Solar energy can be used easily to provide 
this warm temperature source for the heat pump evaporator 
since only a low grade temperature (50-100 degrees fF) is 


required. 


Over the past several years many papers and studies 
have been done on the evaluation of solar energy heat pump 
systems. The earliest research on solar assisted heat pumps 
was begun by the American Gas and Electric Service 
Corporation in 1949 in a5 year study which used a high 
temperature storage system with the heat pump coupled 
directly to the collector[2]. In 1952 R.~ Ce Jordan and S. 
Le Threlkeld published a paper in which the potential of 
this system on the basis of some experiments was 
described[({3].- Since that time the solar assisted heat pump 
system has been improved upon and extensively studied. The 
effects of geographical location and variations in the 
system design or component design have been addresged in the 


present literature [4-6]. 


The study described in this paper will analyze _ the 
performance of a solar assisted water-to-air heat pump 
residential heating system in State College, Pennsylvania 


(latitude 40.5 degrees N). A realistic residence and solar 
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assisted water-to-air heat pump system are modeled to study 
the performance characteristics of such a system in a 
northern climate. The transient simulation computer code 
TRNSYS developed by the University of Wisconsin and as 
modified at the Pennsylvania State University is used. The 
system is studied over a one month winter period, December, 
using actual hourly weather data for State College. The 
system is analyzed for both the cloudiest and clearest 
December weather recorded in the last 30 years. The 
collector area and storage tank capacity are varied and the 
effects on system performance are studied- Other parameters 
varied and studied are the collector/storage tank heat 


exchanger effectiveness and the room temperature. 
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CHAPTER 2 


2-1 Description of Residential Structure 


For investigation and modeling purposes’ the~ solar 
assisted heat pump system was assumed to be installed in a 
single family residence. It is a two story structure with a 
full basement and attached two car garage- The living area 
is approximately 2200 square feet with 8 foot ceilings 
throughout. The solar collectors are located on the south 
slope of the roof and are tilted at 50.5 degrees (latitude + 
10 degrees), the recomended [2] optimum angle for heating 
applications. A floor plan of the residence is shown in 
Figure |]. The residence is located in State College, 
Pennsylvania, a small university oriented city oeituated in 


the center of the state. 


An energy conserving design has been utilized to 
minimize the required heating load. The house walls are 
built of wood frame construction with four inch face brick. 
All windows are double pane glass and ae storm door is 
installed on the front door. Three and one half inches of 
R-1l1 insulation is installed in the walls and six inches of 
R-22 insulation is placed in the second floor ceiling. A 
more complete list of the construction details can be found 


in Appendix A. 
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FIGURE 1 = Floor 
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plan of residence utilizing a 
assisted heat pump system in 
State Collepe, Pa. 
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In order to size and select the appropriate heat pump 
it was necessary to calculate a design summer cooling load 
as well as the winter heating load. The heating load was 
calculated using the general procedures set forth by 
ASHRAE[3] using ASHRAE tabulated data. The resulting design 
heat loss was approximately 32222 Btu/hr. From this’ the 
building UA (thermal transmittance X exposed area) was found 
to be 511 Btu/hr-F which indicates avery well insulated 
house.e- The cooling load was determined using the ASHRAE 
simplified procedures for residential calculations [4] and 
ASHRAE tabulated date. The cooling load requirement was 
found to be approximately 1-4 tons. Summary calculation 
sheets for both the heating and cooling loads can be found 


in Appendix A. 


2-2 System Description 


A schematic diagram of the residential solar assisted 
heat pump heating system is shown in Figure 2. The major 
components of the system are a flat plate solar collector, 
storage tank and a water-to-air heat pump. Other equipment 
includes pumps, heat exchanger, a relief valve and auxiliary 
electric heating which is an integral part of the heat pump. 
The major components are discussed in more detail in the 


following paragraphs. 
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FIGURE 2 = Schematic of solar assisted 
heat pump heating system 
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Ae Flat Plate Collector - The flat plate solar collector 
chosen for the simulation was a Pittsburg Plate Glass(PPG) 
Model CD414. The PPG collector features weather’ sealed 
construction with two tempered glass cover plates and an all 
copper absorber plate with a non-selective coating. The 
collector is especially designed for applications where the 
collector area is limited or where large temperature 
Mieeerentials exist between the ambient air and the absorber 
plate. Each collector panel is approximately 20 square 
feet. The collector array uses a 40% ethylene glycol-water 
solution as the circulated heat transfer medium, This 


provides protection from freezing to -12 degrees F. 
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ADDITIONAL PARAMETERS 


Wind Velocity: 0-5 mph 

Fluid: 40% ethylene glycol solution 
Cp = .792 Btu/lb-=F 

Flowrate = .025 gpm/sqft 

Tube Diameter = .5 in 

Tube Spacing = 3 in 
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FIGURE 3 - Efficiency of the PPG flat plate 
solar collector 


The performance of this flat plate collector can be 


described by a number of parameters- Perhaps the most basic 


performance parameter is the collector efficiency( ) which 


for the PPG collector is 


=.69-.7((T£i-Ta) /HT] (1) 


where Tfisinlet fluid temperature(degrees F) 
Ta =ambient air temperature(degprees F) 

HT incident radiation per unit area on 
the collector surface(Btu/hre-saqft ) 
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The efficiency equation is shown graphically in Figure 
3 along with a number of other parameters describing the 
construction of the PPG flat plate collector. These 
parameters include the transmittance (% of the cover, the 


absorptance(&) of the absorber plate and the emmissivity (€) 


of the absorber plate. 


Three other useful parameters describing the PPG 
collector are the Collector Efficiency Factor(FR), the 
Collector Overall Energy Loss Coefficient (UL) and_ the 
Collector Geometry Efficiency Factor(F’). These parameters 
are related to the Rate of Useful Energy Collected per Unit 
Area of Collector(qu) according to the following expression 
developed by Hottel and Woertz(5] and extended by 


Whillier{6] and Bliss{/7]. 


qu@FR(HT(TK)-UL(T£i-Ta)) (2) 
Since ML =qu/HT equation(2) can be rewritten as 
=FR(TX)-(FR*UL(T£i-Ta)/HT) (3) 


Comparing this to the collector efficiency equation for 


the PPG collector, it follows that 


FR(TX)=.69 and FR*UL@=. 7 (4) 


== 





Substituting known values for the transmittance and 


absorptance 


FRe.69/(€.72)(.95)#1.0 (5) 


thus 


UL=.7 Btu/hr-sqft-F (6) 


Finally, using UL, the Collector Geometry Factor(F 9 
can be obtained from graphs developed by Duffie and 


Beckman[8]). The value of F’ for the PPG collector is .99. 


B. Storage Tank - An tinsulated steel underground tank 
containing water is used for thermal storage. The overall 
thermal transmittance(U) of the tank is .05 Btu/hresqft-F. 
The tank is cylindrical in shape with a constant length to 
diameter ratio of three(as the volume is varied in the 
simulation)- A tube-in-shell heat exchanger with an assumed 
effectiveness of /0Z is used to transfer heat from the 


collector loop to the storage tank. 


C. Water-to-Air Heat Pump -=- An American Air Filter(AAF) 
water-to-air heat pump Model HW19 (1-5 tons) was selected 
for the residencee- As is standard practice the heat pump 
was selected to meet the cooling load requirements'7 of the 
building. It is recognized that in harsh weather the heat 
pump will be unable to provide 100% of the heating load 


therefore the AAF heat pump is equipped with a parallel 
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auxiliary electric resistance heater to provide the deficit. 
The AAF Model HW19 water-to-air heat pump is’ rated to 
provide 600 CFM of heated air to the house. The water flow 
on the water side of the heat pump is 5 gpm- Table 1 is the 
heating performance characteristics of the Model HW19 for 


maintaining the room temperature at 70 degrees F. 


TABLE 1 - Heating performance characteristics 
of AAF Model HW19 heat pump 


EVAP HEAT HEAT 
INLET TEMP ABSORBED REJECTED POWER INPUT COP 
(degrees F) (MBtu/hr ) (MBtu/hr) (kw) (MBtu/hr ) 


40 11.5 18.0 1.80 6.14 2.93 
50 13.3 20.7 2-06 7-03 2-94 
60 15.2 2 Sane 2-30 7-85 2-96 
70 17.1 25.9 a | 8.63 3-00 
80 19.0 28.4 2-76 9.42 3-01 
90 20.9 31.1 2-99 10.20 3-05 
100 22.8 33.7 3-22 10.99 3-07 
110 24.7 36.3 3-45 11.77 3-08 


2-3 Operation of System 


The solar assisted heat pump system is designed to 
supply the complete heating requirements of the residence 
with minimal energy consumption. The required heat is 
provided by one of three modes.~ These modes are direct 


heating from storage, heat pump heating. and auxiliary 


electric resistance heating. 
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The collector loop uses a 40% ethylene glycol-water 
mixture to prevent freezing during low ambient temperature 
conditions. The mixture is pumped through the collector/ 
heat exchanger loop at a rate of .025 gpm/sqft of collector. 
The pump tis controlled by a control mechanism which turns 
the pump off when the temperature of the fluid leaving the 
collector is less than the storage tank temperature. The 
pump resumes operation when the temperature of the fluid in 
the collector reaches the storage tank temperature plus 10 
degrees F. This control sequence prevents the loss of 


thermal energy from the tank by reradiation through the 


collector at night or cloudy days when solar insolation 


decreases to a minimum. 


When the house thermostat calls for heat and _ the 
storage tank temperature is greater than or equal to 110 
degrees F water from the tank is pumped through a heating 
coil in the house air duct while the house air fan 
circulates warm air through the house. This is direct 
heating from the storage tank. If heating is required and 
the storage tank temperature is below 110 degrees F_ the 
direct heating mode is disabled and the control valves 
circulate the water through the heat pump heat exchanger. 
The heat pump now provides heat to the house. On occasion 
the capacity of the heat pump will be unable to’ meet the 
house heating requirement in which case the built-in house 


air stream electric resistance heater of the heat pump will 


ee 








provide the additional required heat. If the storage tank 
temperature falls below 40 degrees F and heating is still 
required then the heat pump is shut down and all the 
required heating is provided by the electric resistance 


heater. 


If the ambient air temperature exceeds 65 degrees F the 
heating requirements decrease to zero and all heating modes 


become inoperative. 


2-4 Computer Model 


The solar assisted heat pump system shown in Figure 2 
was modeled using the Transient Simulation Program (TRNSYS) 
developed at the University of Wisconsin. The program 
allows one to make a transient (time dependent) simulation 
of the overall system by interconnection of individual 
component models. The TRNSYS modeling components’ used are 
shown in Figure 4 which is a flow diagram of the computer 


model. 


The reader component is used to input hourly weather 
data.- The simulation is run for one month of the heating 
season. Only ambient temperature and solar insolation 
(measured on a horizontal surface) data are needed for the 
simulation. Actual weather data obtained from the 


Pennsylvania State University weather station for the months 
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FIGURE 4 - TRNSYS computer model flow diagram (Terms 


are defined in Appendix B) 
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of December 1960 and December 1972 are used. For the past 
thirty years of available weather data, these years 
represent the highest and the lowest values for December 


sOlar insolation, respectively, for State College, Pa. 


The radiation processor component converts the solar 
insolation or radiation ona horizontal surface (HH), to 
solar radiation on a tilted surface (HT)- The collector is 
tilted at an angle of 50.5 degrees. By using Mode 3 of the 
radiation processor component the total radiation on a 
tilted surface is a combination of the beam radiation, 
diffuse radiation and reflected radiation. The radiation 
processor calculates these values from the horizontal 
surface radiation, the collector orientation and the ground 


reflectance which is -2 for bare ground. 


The liquid collector-storage module is a single 
cOmponent used to model the solar collector, ethylene 
glycol/water heat exchanger, storage tank and associated 
pumps and relief valve. A schematic of the module is shown 
in Figure 5. The combining of these components’ into one 
subsystem reduces computation time and expense. The 
collector fluid, 40% ethylene glycol solution, is pumped at 
the collector manufacturers recommended rate af .025 
gpm/sqft of collector. Across the heat exchanger’ the 


capacity rate ratio is held at 1.0 [i-e. (mCp)eg = (nCp)w ]). 
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The initial storage tank temperature at the beginning of 
December is assumed to be 70 degrees F. It turns qut that 
any initial temperature in this range gives nearly identical 
results over the one month period. The relief valve is set 


to maintain the tank temperature below 212 degrees F. 


RELIEF = 
VALVE 
STORAGE 
TANK 


FIGURE 5 = TRNSYS liquid collector/storage module 
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The water-to-air heat pump component uses actual heat 
pump performance characteristics. The model is quasi-steady 
state in nature- That is, at any instant in the "transient" 
simulation, the heat pump performs as,if it were at steady 


state, having for the moment constant source and sink 


temperatures”~ and hence a unique heat addition, heat 
rejection and work input. The heat pump is modeled _ to 
Operate in the heating mode only. A parallel auxiliary 
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electric heating element is included to add just enough 

energy to meet the heating load when the heat pump capacity 

is less than the load or when the’ storage tank temperature 
\ 

falls below 40 degrees F. Figure 6 is a schematic of the 


water-to-air heat pump component. 
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FIGURE 6 = TRNSYS water-to-air heat pump component 


The heating load component determines the heating load 
of the residence by use of the energy/degree day concept. 
The concept is extended to estimate the ENE Ole heating 
load by multiplying the heating requirement of the residence 


(energy/degree-hour) by the difference between the required 


room temperature and the outdoor ambient temperature. 


Beg 








The quantity integrator simply keeps a continuous total 
of various system outputs throughout the simulation. It 
works similiar to a KWH meter which continuously totals the 


amount of electrical energy consumed. 


The printer components control the display of the 
system outputs. In this simulation some of the outputs are 


printed at hourly intervals and some at daily intervals. 


Each of the TRNSYS components’ requires ae set of 
defining parameters and inputs to produce an output. The 


parameters describe the operating characteristics of the 


system component being modeled. Many of these parameters 
have been discussed previously. A complete list and 
description of the typical parameters used in this 


simulation for each component can be found in Appendix B. 
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CHAPTER 3 


3.1 Presentation and Discussion of Results 


The collector efficiency, which is a measure of the 
useful solar energy being transferred to the heating systen, 
is shown in Figures 7 and 8 for various values of collector 
area and storage capacity. The collector efficiency is 
maximum where the collector area is small and the storage 
capacity is large. This can be shown also by the collector 
efficiency equation when one realizes that the inlet fluid 
temperature to the collector (from storage) will be a 
minimum when the ratio of collector area to storage capacity 


is a minimum. In other words looking at equation (1) again, 


y = .69-.7((Tfi-Ta)/HT], (1) 


low values of Tfi (resulting from a small collector grea to 
storage capacity ratio) will reduce the second term on the 
right hand side resulting in a larger value for the 
collector efficiency. The collector efficiency increases 
for increases in the storage capacity and decreases in the 
collector area since both of these trends decrease the 
temperature of the fluid entering the collector. Figure 7 
also shows that increasing the storage capacity beyond about 
2 gal/sqft of collector area is of little benefit since the 


collector efficiency remains relatively constant. 
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The effect of weather on the collector efficiency can 
be seen by comparing the curves for 1960 and 1972 in Figures 
7 and 8. The 1972, or minimum solar insolation, curves 
reflect higher collector efficiencies for identical system 
parameters (collector area and storage capacity). This 
result is not as easy to predict from equation (1) since all 
three variables on the right Side change. From the 
graphical results, however, it can be concluded that for 
poor weather the predominant effect on the efficiency 
equation is a decrease of the second term on the right hand 
side. This implies that the magnitude of the drop in 
entering fluid temperature must be larger than’ the drop in 
solar insolation, indicating relatively low entering fluid 
temperatures (storage tank temperatures) would be expected 
for poor weather conditions. Figure 7 also shows” that 
weather extremes have a negligible effect on the collector 
efficiency of systems with small collector area and a 
Significant effect on the collector efficiency of systems 


with large collector area. 


Values of the System Coefficient of Performance (SCOP) 
for various collector area and storage capacity 
configurations are shown in Figures 9 and 10. The SCOP is 
defined by the following equation: 


total heating load 
SCOP @ wwm nnn mene eee mt mm ew eee eee ee eee  ateadaadeadeentied (7) 


electric energy used by heat pump & aux heater 
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Following conventional practice, the energy required for 
fans and pumps is considered negligible in comparigon with 
other energy inputs. Since the total heating load 
requirement for the residence is independent of the heating 
system configuration, the SCOP will depend only upon the 
amount of energy used to run the heat pump and the energy 
required for the auxiliary heater. The required energy 
input from these two sources is reduced, however, when there 


is a contribution from direct heating. 


In Figures 9 and 10 it can be seen that all the curves 
approach a SCOP of 3-0 as the overall system size (collector 
area and storage capacity) increases. This is because as 
the storage capacity to collector area ratio increases the 
system becomes more and more dependent upon the heat pump 
alone (i-e. the direct heating contribution and the 
auxiliary heat contribution to the heating load decrease). 
Specifically, the curves converge to a value of 3-0 since 
that is the average heating COP of the AAF heat’ pump over 


its range of evaporator inlet (storage tank) temperatures. 


Figures 9 and 10 both show that the maximum SCOP is 
obtained with a large collector area and a_e small storage 
capacity. The factor responsible for driving the SCOP past 
3-0 (the average heating COP of the heat pump) is the direct 


heating contribution. Without a direct heating contribution 
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the SCOP will be limited to a maximum value equal to the 
heating COP of the heat pump. Also from Figure 1Q@ it is 
seen that for a SCOP greater than 3.0 the’ solar assisted 


heat pump system must have as a minimum a collector area 


greater than about 650 sqft. 


The collector area of 650 sqft in Figure 10 is also 
important in that it is the approximate collector area for 
which all the different storage capacity curves resulting 
from 1960 weather data intersect each other at a SCOP equal 
to 3.0. The significance of this intersection is that it 
represents the point at which the total residence heat 
requirement is provided by the heat pump and system 
performance is independent of the storage capacity. This 
means that at 650 sqft of collector area the energy input 
from the solar collector array is equal to the energy 
withdrawn from the storage tank by the heat pump, resulting 
in neither the benefit of a direct heating contribution or 
the requirement for auxiliary heat. By extrapolation it is 
estimated that a similiar intersection point will occur for 
the storage capacity curves resulting from 1972 weather data 
at about 1800 sqft. Since 1960 and 1972 represent’ the 
extremes of solar insolation it is logical to expect that 
for the average year. the intersection point will be 
somewhere between 650 and 1800 sqft, say 1200 sqfte Thus, 


based on a year of average weather the solar assisted heat 
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pump system for State College would require a collector area 
of 1200 sqft or greater to obtain a SCOP greater’ than 3.0 


(the heat pump heating COP). 


For collector areas below the intersection point the 
residence heat load is supplied essentially by the heat pump 
and auxiliary heat and the SCOP increases as the storage 
capacity increases. For collector areas above the 
intersection point the residence heat load is supplied by 
the heat pump, auxiliary heat and direct heat and the SCOP 


decreases as the storage capacity increases. 


This situation suggests that a two tank storage system 
would be beneficial.- A smaller high temperature tank would 
be charged first by the solar collectors and be used for 
direct heating. A second larger tank in the medium 
temperature range would be charged with the excess’ solar 
energy and be used as a heat source for the heat pump. This 
combination would allow the direct heating contribution to 
be maximized. Whether or not the advantages of a two tank 
system outweigh the added cost and complexity would have to 
be studied in more detail.- Figures 9 and 10 also indicate 
that for a single storage tank system the SCOP is’ more 
strongly influenced by the collector area than by the 
storage capacity, since in general the SCOP increases more 
for an increase in collector area than for an increase in 


storage capacity. 
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Comparing the curves resulting from 1960 weather data 
(good) and 1972 weather data (poor) in Figures 9 and 10 
shows that the weather data used for a_ solar assisted heat 
pump system simulation has a major effect on the resulting 
system performance characteristics. The SCOP obtained for a 
given collector area~-storage capacity configuration can vary 
by more than 250%. If non-representative weather data are 
used in the design calculations, a system design could 
result which is inadequate for one situation or overdesigned 


for another. 


Figures 11, 12 and 13 show the precentage contribution 
made by the heat pump, direct heat and auxiliary heat to the 
total heating load for a range of collector areas and 
storage capacities. For a given collector area the heat 
pump contribution increases as the storage capacity 
increases. Accordingly, the percentage of the total heat 
load supplied by auxiliary heat decreases. Looking at the 
bar graphs for 1960 and 1972 (representing weather extremes) 
it is seen that poor weather, as intuitively expected, 
increases the auxiliary heat requirements. The amount of 
increase is most severe around collector areas of 600 sqft 
where it averages about 40%. The percentage contribution 
made by the heat pump, for a fixed storage capacity to 
collector area ratio, peaks in the middle range of collector 


areas. As predicted previously by Figures 9 and 10, the 
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heat pump contribution is near 100% for a collector area 
near 650 sqft (for 1960 weather data). At large collector 
areas (1000 sqft, Figure 13) the heat pump contribution is 
smaller for good weather than for poor weather because the 
direct heat contribution lessens the heat requirement to be 


provided by the heat pump. 


Direct heating becomes a significant contribution when 
the collector area is large and the storage capacity is 
small because the large amount of collected solar energy 
only has to heat a small amount of water, resulting in a 
high temperature level in the storage tank. Decreasing the 
collector area causes the percentage contribution of direct 
heating to decrease rapidly so that for small collector 
areas the direct heating contribution is non-existent. For 
a sufficiently large collector area, increasing the storage 
capacity will decrease the direct heating contribution. 
Also, poor weather will decrease the likelihood of any 
direct heating contribution. In Figure 13, for a storage 
capacity of 1 gal/sqft of collector, the direct heat 
contribution for poor weather (1972) is only 52% compared to 
46% for good weather (1960). As indicated by previous 
figures, the bar graphs of Figures 11, 12 and 13 show that 
the best performimg configuration of the solar assisted heat 
pump system is a large collector area with a small storage 
capacity, mainly because of the advantage of more direct 


heating. 
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Table 2 lists the energy usage of a system with 1000 
sqft of collector area and a storage capacity of 1 gal/sqft 
of collector. Energy data is given for both "gaod" and 
"poor" weather. It is clear from the percent deviations 
listed that weather has a major impact on the system 
performance and energy requirements. Note that when the 
direct heating contribution falls off during poor weather 
the majority of the difference is made up by auxiliary heat. 
Several comparisons can be made from Table 2 besides the 
deviations shown. For instance, assuming an electric rate 
of $.035/kwhr for State College and the system configuration 
given in Table 2, it would cost between $28 and _ $73 
(depending upon the weather) to heat the house for the month 


of December. 


TABLE 2 - A one month comparison of system performance 
characteristics for a solar assisted heat pump 
System in State College, Pa. ; Collector area = 
1000 sqft and water storage capacity = 1000 gal 


Dec 1960 Dec 1972 Percent 


Weather Weather Deviation 

Data Data From 1960 
Collector Efficiency (Z) 36 44 MILE 
Sys Coefficient of Performance 5-14 1.90 -~63.0 
Total Energy Collected (kwhr) 3648 1778 -51.2 
Total Heating Load (kwhr) 4161 3950 -5.1 
Heat Pump Contribution (kwhr) 2151 2508 16.6 
Direct Heat Contribution (kwhr) 1902 193 -89.9 
Aux Heat Contribution (kwhr) 100 1242 1142 
Energy Input to Heat Pump (kwhr) 709 84 1 18.6 
Energy Input to Aux Heat (kwhr) 100 1242 1142 
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Since the collector array is mounted on the roof of the 
residence, in reality an upper limit of collector array 
size is set by the area of the roof. Taking this into 
consideration a collector area of 600 sqft was used for the 
parametric studies of room thermostat setting and ethylene 
glycol/water heat exchanger size. Figure 14 shows for a 
collector area of 600 sqft the effect on SCOP of lowering 
the room temperature to 65 degrees F as the storage capacity 
is varied. It is interesting to note that the SCOP is 
increased when the room temperature is maintained at 65 
degrees F versus 70 degrees Fe. Recalling that the SCOP is 


defined by equation (7) as 


total heating load 


electric energy used by heat pump & aux heater 


one might expect that since the heating load tis reduced by 
about 13% the SCOP would decrease, but the reduction in 
energy input to the system is even greater. The heat pump 
input energy is reduced by 13% and the auxiliary heat input 
is reduced to nearly zero resulting in a net increase in the 
SCOP. As shown by the nearly horizontal curves of Figure 
14, and as previously stated for a collector area of 600 
sqft, the storage capacity size has little effect on the 


SCOP. 


Table 3 shows how lowering the room temperature to 65 
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TABLE 3 = Comparison of energy requirements for room 
temperatures of 65 F and 70 F for a solar assisted 
heat pump system in State College, Pa. 


DEC 1960 WEATHER DATA 


Collector Area = 600 sqft 
Storage Capacity = 6000 gal 


Room Temp Room Temp 

65 F 70 F 
Total Heating Load (kwhr) 3604 4161 
Heat Pump Contribution (kwhr) 3598 4120 
Direct Heat Contribution (kwhr) 0 0 
Aux Heat Contribution (kwhr) ~03 34 
Energy Input to Heat Pump (kwhr) 1193 1374 
Energy Input to Aux Heat (kwhr) ~ 03 34 


DEC 1972 WEATHER DATA 


Collector Area = 600 sqft 
Storage Capacity = 6000 gal 


Room Temp Room Temp 

65 F 70 F 
Total Heating Load (kwhr) 3393 3950 
Heat Pump Contribution (kwhr) 2622 2683 
Direct Heat Contribution (kwhr) 0 0 
Aux Heat Contribution (kwhr) 765 1263 
Energy Input to Heat Pump (kwhr) 888 909 
Energy Input to Aux Heat (kwhr) 765 1263 


one a 





degrees F affects the energy inputs of the system. A 


storage capacity of 6000 gallons was chosen since Figure 10 
shows that for 600 sqft of collector area the maximum SCOP 
is obtained with a storage capacity of 10 gal/sqft of 


collector area. The economics of lowering the room 


thermostat setting can be obtained from the values in Table 


3. Again assuming an electric rate of $.035/kwhr for State 


College it is seen that for the system configuration given 
in Table 3 approximately $8 to $18 can be saved (depending 


upon the weather) for the month of December by lowering the 


thermostat setting 5 degrees. 


Tables 4 and 5 are the system performance 
characteristics for good (1960) and poor (1972) weather 
resulting from variations in the ethylene glycol/water heat 
exchanger effectiveness. If the capacity rate ratio is held 
constant at a value of 1.0 (ieee. Cmin/Cmax = 1.0) the heat 


exchanger effectiveness (E) is defined by the following 


equation: 


E = NTU/(1+NTU) (8) 


where NTU = Number of transfer units 
(size parameter) 


Thus by holding the capacity rate ratio constant at a value 
of 1.0, variations in the heat exchanger effectiveness (E) 


actually represent variations in the size of the heat 
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TABLE 4 - 
pump system 


effectiveness based 
Pa. 


State College, 


STORAGE COLLECTOR 
TANK 
TEMP(F) Ce) 


HXGR 
EFF 


(Collector Area=600sqft 


System performance factors of 
for varying values of 


EFFICIENCY SYSTEM 
COP 


solar assisted heat 
heat exchanger 


on DEC 1960 weather data for 


CONTRIBUTION TO HEAT LOAD (2) 
HEAT DIRECT AUX 
PUMP HEAT HEAT 


Storage Capacity#=600gal) 


- 6 115.8 44.3 2-68 73 14 13 

-7 117.3 44.7 Ziel 72 15 13 

8 118.3 45.0 2-74 fal 16 13 
(Collector Area=600sqft Storage Capacity#=#2400gal) 

° 6 87.1 47.7 2-93 99 0 l 

27 88.8 48.0 2-94 99 0 l 

8 90.0 48.3 2-95 99 0 1 


(Collector Area600sqft 


Storage Capacity#=6000gal) 


«6 76.7 48.5 2.95 a9 0 l 
2/7 78.0 49.0 2.96 99 0 l 
8 78.9 49.4 2-96 99 0 ] 


TABLE 5 = System performance factors of 
for varying values of 


pump system 


effectiveness based 
Pa. 


State College, 


STORAGE COLLECTOR 
TANK 
TEMP (F) 


HXGR 
EEF 


(4%) 


(Collector Area=600sqft 


EFFICIENCY SYSTEM 
COP 


solar assisted heat 
heat exchanger 


on DEC 1972 weather data for 


CONTRIBUTION TO HEAT LOAD (2) 
HEAT DIRECT AU X 
PUMP HEAT HEAT 


Storage Capacity#=600gal ) 


- 6 40.2 48.8 1.47 46 l 53 

27 39.8 49.4 1.48 46 ] 53 

8 40.1 49.8 1.48 47 l 52 
(Collector Area=600sqft Storage Capacity#=2400gal) 

- 6 40.0 52.6 1.62 57 0 43 

| 40.0 53.3 1.63 58 0 42 

8 40.0 53.8 1.63 58 0 42 
(Collector Area600sqft Storage Capacity#=6000gal) 

- 6 40.0 52-0 1.81 67 0 33 

77 40.0 52.7 1.82 68 0 32 

8 40.0 53.3 le 68 0 a2 
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exchanger.e. The smaller the value of E the smaller is the 
size of the heat exchangere AS can be seen from Tables 4 
and 5 varying the heat exchanger effectiveness (i.e. size) 
of the ethylene glycol/water heat exchanger has negligible 
effect on the performance of the system. The SCOP increases 
very slightly as the heat exchanger size increases 


regardless of the weather for the effectiveness ' range 


studied. 
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CHAPTER 4 


ie 1 Conclusions 


The maximum collector efficiencies are obtained for 
system configurations with small collector areas and large 
storage capacities because the storage tank temperature and 
thus the temperature of the fluid entering the collector is 
kept low. Similarly, the collector efficiency of a given 
system is improved during periods of poor’ solar insolation 
due to the reduced entering fluid temperature. The 
collector efficiency increases as the storage capacity 
increases and as the collector area decreases since the 
combination of these result in a lower temperature of the 
fluid entering the collector. Weather extremes have little 
affect on the collector efficiency of systems with small 
collector area, however, the affect on the collector 
efficiency is significant for systems with large collector 
area. The study results also show that for the solar 
assisted heat pump system a storage capacity beyond 2 
gal/sqft of collector area is of little benefit since the 


collector efficiency remains relatively constant. 


For a solar assisted heat pump system the SCOP is 


limited to the heating COP of the heat pump unless’ the 


combination of good solar insolation and adequate collector 
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area result in a heating contribution from direct heating. 
The SCOP will be largest for systems with large collector 
areas and small storage capacities. There is a unique 
collector area at which the SCOP of the solar assisted heat 
pump system will equal the heating COP of the heat pump 
(i-e. 100Z% of the total residence heating requirement is 
provided by the heat pump) for a wide range of storage 
capacities. This unique collector area is dependent upon 
the solar insolation received and will vary in value 
according to the weather. For the solar assisted heat pump 
system modeled for State College in this study it is 
estimated that about 1200 sqft of collector area would be 
required to provide 100Z heat pump heating under’ average 


December weather conditions. 


The SCOP of the solar assisted heat pump system 
increases with increasing collector area and in general is 
more strongly influenced by the collector area than by the 
storage capacity. For a system with a collector area 
greater than that resulting in 100% heat pump heating, the 
residence heating load is supplied by contributions from the 
heat pump, auxiliary heat and direct heat and the SCOP 
decreases as the storage capacity increases. For a system 
with a collector area smaller than that resulting in 100% 
heat pump heating, the residence heating load is supplied 


essentially by the heat pump and auxilary heater and the 
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SCOP increases as the storage capacity increases. This 
situation suggests the use of a two tank storage system with 
a high temperature tank being used to maximize the direct 


heat contribution and a medium temperature tank for use as 


the heat pump heat source. 


For a given collector area the percentage contribution 
of the heat pump increases as the storage capacity increases 
whereas the percentage contribution of the auxiliary heater 
decreases. Direct heating becomes a significant 
contribution to the total heating load for solar assisted 
heat pump systems with large collector areas and small 
storage capacities since the storage is heated to a high 
temperature level. Decreasing the collector area causes the 
percentage contribution of dtrect heating to decrease 
rapidly so that for small collector areas the direct heating 


contribution is non-existent. 


When designing or simulating solar assisted heat pump 
systems it is essential that accurate and realistic weather 
data be used. The weather data used has a significant 
effect on the performance characteristics of the system. 
“For the system modeled in this study the SCOP varied by as 
much as 250% for some collector/storage configurations when 
"poor" versus "good" weather data were used. The percentage 


contributions of the heat pump, auxiliary heater and direct 
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heating to the total heating load are similarly affected as 
was shown in Figure 13. A simple economic analysis of a 
system in State College with 1000 sqft of collector and 1000 
gallons of storage shows that monthly heating costs can vary 


from $28 to $73 depending upon the weather. 


Federal energy conservation guidelines call for 
reducing thermostat settings to 65 degrees F for winter 
heating. This study found that reducing the room 
temperature from 70 degrees F to 65 degrees F increased the 
SCOP approximately 5% and reduced significantly the amount 
of auxiliary heat prequired- For a system configuration of 
600 sqft of collector area and 6000 gallons of = storage 
capacity a simple cost analysis shows that in State College 
approximately $8 to $18 per month (depending upon the 


weather) can be saved by lowering the thermostat setting 5 


degrees. 


The performance of the solar assisted heat pump system 
is not significantly affected by the effectiveness (i.e. 
size) of the ethylene glycol/ water heat exchanger. Since 
this is the case, it appears that the heat exchanger could 


be of relatively simple design and low cost. 


When reviewing the results presented in this study it 


must be remembered that several assumptions have been made 
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in modeling and simulating the operation of the residential 
solar assisted heat pump system. The majority of these 
assumptions are inherent in the TRNSYS components but others 
are externally controlled. For instance the effect of wind 
velocity has been assumed negligible which in a "real world" 
installation could reduce the solar energy collected by the 
collectore Also a ground reflectance value for bare ground 
has been assumed. In a "real world" installation’ the 
reflectance value might be higher if the ground was snow 
covered resulting in an increase of reflected solar 
radiation received by the collector. Overall, however, it 
is expected that the model performance will be somewhat more 


optimistic than that of a "real world" installation. 


Although domestic hot water heating was not included in 
the solar assisted heat pump model of the _ residential 
heating system, the collector array sized for the solar 
assisted heat pump system would most likely provide all of 
the domestic hot water requirements in the summer months. 
During the heating season the system would be able _ to 
provide preheating of the domestic hot water much of the 


time. 


At first glance the TRNSYS modeling program seems very 


complicated but actually the building block modeling 
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technique is fairly straightforward. Individual system 
components (i.e. pumps, solar collectors, storage tank, 
etc.) are described by component models which are connected 
together to form a complete system. It is better to keep 
the model as simple as possible since problems, such _ as 
failure of the differential equation solution techniques to 
converge within the component models, are likely to result 
if the model is too elaborate.- Some modes of the TRNSYS 
components use empirical equations which have been deyeloped 
in SI units requiring the entire simulation to be run in SI 
units if these component modes are used. If these component 
modes can be avoided in the model the more familiar English 
units can be used. The TRNSYS simulations require hourly 
weather data as input-e The assembly and handling of such 
data can become quite cumbersome if punched cards are being 
used, especially if the simulation is run for a full year. 
This problem can be alleviated by using a tape input. The 
computer time required for running a full year simulation of 


a sOlar assisted heat pump system (based on this one month 


run) will be about 200 seconds and cost approximately $60 if 


output is printed for each hour. 


4.2 Recommendations 


To obtain a complete picture of the solar assisted heat 


pump system, a detailed economic analysis of the system 
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alternatives should be performed. This would determine the 
most cost effective combination of such things as collector 
area and storage capacity. An economic analysis would also 
be an effective method to compare the solar assisted heat 


pump system with conventional heating systems. 


An entire year or at least a complete heating season of 
weather data should be used when simulating the solar 
assisted heat pump system rather than a single month. This 
will give a more realistic picture of the system’s seasonal 
performance characteristics. Limiting the weather data to 
the heating season will significantly decrease the amount of 
data and reduce computer costs.- It is important however 
that the weather data be carefully selected to provide a 


realistic representation of the location being studied. 


Future studies on the solar assisted heat pump system 
should explore the feasibility and operating characteristics 
of a two tank storage system for optimizing the output of 
the solar assisted heat pump system. Also the results 
presented in this report and any future simulations using 
the TRNSYS program should be compared with other modeling 
techniques such as those being developed by the 
Architectural Engineering Department of the Pennsylvania 


State University. 
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APPENDIX A 
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BUILDING CRITERIA 


Indoor Design Temperature = 70F 


Outdoor Design Temperature = 7°F 


CONSTRUCTION DATA 


Wall: 4"' Face Brick 
3/4'' Air Gap (Mean Temp. = 50, Temp. Diff, = 10, E = . 82) 
Building Membrane (Permeable Felt) 
1/2'' Plywood (Douglas Fir) 
3 1/2"' Fiberglass Insulation 
1/2'' Gypsum Wall Board 
2x 4 Stud Wall Construction, 16''0O,C, 


Glass: Area is given on plan; e.g. A = 15 
All Double Pane 
(1/4'' Air Gap, Wood Sash, 80% Glass) 


Doors: Sis aay Ge 
1 1/2'' Solid Wood 
Metal Storm Door (Front Door only) 


Ceiling: 6'' Fiberglass Insulation 
1/2"' Gypsum Board 
2x 6 Ceiling Joists, 16'' O.C, 


ass ee 





U-Values: 1977 Handbook of 




















Walls: R Fundamentals Reference 
fice OF Table No. 
Outside Air Film (15 mph wind) - yay) [-A 
4'' Face Brick Zag is 32-A 
3/4'' Air Gap O 2Z AZ 7. 
(Mean Temp. = 50, Temp. Diff.= 10, 
=o) 
Building Membrane (Permeable Felt) 0b 22, [a 2- AA 
1/2"' Plywood (Douglas Fir) .6O2 22,13 2-A 
3 1/2"' Fiberglass Insulation .00 2 ele 3-A 
1/2'' Gypsum Wall Board as 72 
Inside Air Film LS ya. ty 
TOTAL RESISTANCE [4,43 
FINAL U-VALUE -O6F 
Glass: U 
Double Pane 
(1/4"' Air Gap) eS 22.24 8-A 
Adjustment for Wood Sash, 80% Glass 
(Ux adjustment) ; IS 22.24 S-C 
FINAL'U-VALUE 55 
Door: U 
1 1/2" Solid Wood .47 
with metal storm door ,33 22.25 | 
FINAL U-VALUE ,33 
Ceiling: R 
Air Film (Still Air, E = .90) .of raCoea hI /-A 
6"' Fiberglass Insulation 2,00 22.14 2. 
1/2'' Gypsum Wall Board YP os 3- 
Inside Air Film b/ ie 7 = 
TOTAL RESISTANCE 23.67 _ 
FINAL U-VALUE Fz 
Kitchen Wall Against Garage: R 
Air Film (Still Air, E = .90) , 63 / [-A 
3 1/2" Fiberglass Insulation £ f : 
1/2'' Gypsum Wall Board Ss rae | <- 
Inside Air Film 6¥ zzZ.1/ ~ 
TOTAL RESISTANCE [2.8 
FINAL U-VALUE .O7 
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HEAT LOSS SUMMARY 


Location (City) 
Design IDBT 2 70F 


Design ODBT = ‘his F 

















Type of Building Z S70KY RESIDENTIAL 
STRUCTURAL AREA U-VALUE UA 
COMPONENT SQ FT) 
S7 [io0R =_ ee 
WALLS TAZ 069 65,00 
GLASS a SS JULES 
FLOOR ASSUME We HEAT LOSS 
_ Doors 8) i tis 
INTERIOR WALL [03 - 
ECON LOOk 
WALLS PAGE 069 77.42 
_ Se er 4/05 
CEILING [os 042 46.4 
Sub TOTAL 
17307 FT? 
INFILTRATION: No. of Air Changes /44 Gr 112 Jue 
Qs 1,1CCFm )AT 
POSITIVE VENTILATION: 0 CFM 
TOTAL 


HEAT LOSS 
10 = ODB1I 


Design Temp. Diff. = 


= 32222 
70-7 


=5a5 


STATE COUuleGe fy 


"KF ind-out temp ditt 














Btuh 


63°F. 


HEAT LOSS, 
Btuh 


_ 405 

4886 

R/S Se 
520 Ate 


= 15 // 





COOLING LOAD CALCULATION 


SOURCES: Chap. 22, 1972 ASHRAE Handbook of Fundamentals 
Chap. 25, 1977 ASHRAE Handbook of Fundamentals . 


Location (City) STATE: COLLEGE, fp 


Design Conditions 75 IDBT, 87 ODBT 
Daily Range 250 E 


Design Equivalent Temperature Differences: 


Brick Walls Soe Ceiling & Roofs 36 F 

Wood Doors /S.6 EF Floors O F 
Design Window Factors: . 

North /3 Stub/ft“ West 43_Btuh/ft® 

East 43 Beuh/ft? South 2/ Btuh/fe® 





= 


COOLING LOAD SUMMARY 


Walls ESE Btu 
Windows 24 lb - 
Doors 32/ o 
Ceiling Loe fal . 
Floor O . 
Occupants 70O 
Cooking . [Zoo s 
ALL ALen) (. 7 brut [er >) Infiltration Ai /4 ‘ 
Total [2304 Btuh, Total Sensible Heat (TS 


7 ee _ ae a 


Total Heat = 1.3 (TSH) = //7Z297S Btuh 


Equipment Size = ee = de 4 4- tons 


aoe 
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DEFINITION OF TERMS USED IN TRNSYS COMPUTER MODEL FLOW DIAGRAM 


DELTAU 


H 

Hb 
Hd 
HT 
HB 
HD 


Hl 


H2 


H3 


IGAM 
Mi 


Mo 


MODE 
OPHRS 


QA 1 


QABS 


QAC 


QAUXH 
QDH 
QDUMP 


QENV 


change in storage tank internal energy 


radiation (total) on a horizontal surface 


radiation (beam) on a horizontal surface 


radiation (diffuse) on a horizontal gurface 


radiation (total) on a tilted surface 


radiation (beam) on a tilted surface 


Solar radiation (diffuse) on a tilted surface 


Fraction of timestep that heat pump operated in 


water-air heating mode 


Fraction of timestep that heat pump operated in 


air-air heating mode 


heat pump operated in 


Fraction of timestep that 
cooling mode 


Heat pump master control integer; Q#off, l=on 
Flowrate of water entering heat pump evaporator 


Flowrate of water returning from heat pump 
evaporator 


Mode of heat pump operation 
Total time heat pump operated in heating mode 


Instantaneous energy absorbed from liquid source in 
water-air heating mode 


Total energy absorbed by heat pump 


Instantaneous energy removed from load in cooling 
mode 


Auxiliary energy delivered to room in heating mode 
Energy delivered to room in direct heating mode 
Rate at which energy is discarded by relief valve 


Rate of energy loss from storage tank to 
Surroundings 
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QL 
QLOAD 


QR1 


QR 2 


QRC 


QREJ 


QSHRTC 


QTANK 


TENV 


TOTQL 


TOTQU 


TOTRAD 


TSTOR 


WAC 


W AH 


WORKIN 


AE 


Instantaneous room heating load 


Hourly heating load 


Instantaneous energy delivered to load in water-air 
heating mode 


Instantaneous energy delivered to load in agir-air 
heating mode 


Energy rejected to sink in cooling mode 
Total energy delivered by heat pump 


Cooling capacity required but not met by heat pump 


Rate of energy delivery from storage tank to heat 
pump 

Rate at which energy is transferred to storage tank 
Hourly useful energy collected 
Total work input to heat pump 


Hourly incident solar radiation on collector 


Temperature of storage tank 

Ambient temperature 

Temperature of ground surrounding storage tank 
Temperature of water entering heat pump evaporator 


Temperature of water leaving heat pump evapprator 


Total heating load 

Total useful energy collected 

Total incident solar radiation on collector 
Hourly storage tank temperature 


Instantaneous work input to heat pump in cooling 
mode 


Instantaneous work input to heat pump in heating 
mode 


Total work input to heat pump in heating mode 


Change in internal energy of water in storage tank 


267 = 





ae 


a 





— 
- Angle of incidence of solar radiation on a tilted 


! ‘e urface 


25.2 








TYPICAL TRNSYS COMPONENT PARAMETERS 


(aA) UNIT 1, 


TYPE 9 


Parameters 


] 
2 


3 


4 
> 


N on 
At - 


i - 


Mi - 
Ai - 


READER 


Description 


number of values to be read 
time interval at which data is 
provided 

position of value on card for 
which units are to be converted 


multiplication factor for that value 3.686 


addition factor for that value 


0 


(parameters 3,4&5 convert Langleys to Btu/hr-sqft) 


6 


(B) UNIT 2, 


7 - 


logical unit number of data file 


TYPE 16 RADIATION PROCESSOR, MODE 3 


Parameter 


sou & WP = 


(Cc) UNIT 3, 


Parameter 


] 
2 
3 
4 
5 


On 


& 


Description Value 
mode 3 
n =- day of year at start of simulation 335 
g - latitude 40.5 deg 
Ss - slope of collector surface 50.5 deg 
Y - orientation angle of collector 0 deg 
Sc - solar constant 429 Btu/hr-sqft 
P - ground reflectance 22 
TYPE 21 LIQUID COLLECTOR=-STORAGE MODULE 
Description Value 
A - collector area 600 sqft 
Mc = collector fluid flowrate 7500 lb/hr 
Cpc- collector fluid specific heat 792 Btu/lb-F 
Mh = heat exchanger to tank flowrate 5940 lb/hr 
Cph- specific heat of fluid in 1.0 Btu/1lb-F 
storage tank 
F’ = collector efficiency factor ~99 
TX - transmittance-absorptance ~684 
product 
UL = collector loss coefficient -7 Btu/hr-sqft-F 
Tmax=- temperature at which tank 212 deg F 
relief valve opens 
€ - collector-tank heat exchanger e 7 
effectiveness 
V =- volume of storage tank 80.1 cuft 
P - tank fluid density 62-4 lb/cuft 
U - tank energy loss coefficient -05 Btu/hr-sqft-F 
r - ratio of tank height to tank 3 


diameter 


=5.0a 





15 To - 
(D) UNIT 4, TYPE 20 WATER-TO-AIR HEAT PUMP 
Parameter Description 
1 CPl -|+ specific heat of ligquid 
entering heat exchanger A 
2 MDOTI1 = mass flow rate of liquid 
entering heat exchanger A 
3 TROOMH = constant room temperature 
for heating mode 
4 TROOMC = constant room temperature 
for cooling mode 
b) TMIN] = minimum liquid source temp 
for water-to-air heat pump 
6 TMIN2 = minimum air source temp for 
air-to-air heat pump 
7 NDATAH = number of equally spaced 
heating data points 
8 NDATAC = number of equally spaced 
cooling data points 
9 ISAVEH |= logical unit number of 
heating data file 
10 ISAVEC = logical unit number of 
cooling data file 
11 TCOOL = minimum ambient temp when 
cooling is allowed 
12 CAIR = MCp product for room air 
flow thru heat exchangers 
B,C and D 
13 EFF - effectiveness of heat 
exchanger B 
14 TSET - minimum liquid source temp 
for direct heating from 
liquid source 
15 ICOOL =- cooling condenser selection 
integer, O=heat exchanger A 
16 THEAT = maximum ambient temp when 
heating is allowed 
(E) UNIT 5, TYPE 12 HEATING LOAD, MODE 3 


Parameter 


l 
2 


3 


initial tank fluid temperature 


Description 


mode 


UA - heating requirement of house 


TR = 


(energy/degree hour) 


room temperature 


-60- 


70 deg F 


Value 
1.0 Btu/lb-F 
2500 lb/hr 
70 deg F 
77 deg F 
40 deg F 
1000 deg F 
8 
8 
31 
32 
1000 deg F 


648 Btu/hr-F 


7 

110 deg F 

0 

65 deg F 
Value 

3 


-511 Btu/hr-F 


70 deg F 





(F) UNIT 6, TYPE 24 QUANTITY INTEGRATOR 


No parameters required 


(G) UNIT 7, TYPE 25 PRINTER (HOURLY) 


Parameter Description Value 
] Tp -~- time interval at which printing l hr 


is to occur 


(H) UNIT 8, TYPE 25 PRINTER (24 hour) 


Parameter Description Value 
] Tp - time interval at which printing 24 hr 


is to occur 


6 1 
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TRNSYS - A TRANSIFNT SEMULATION PROGRAM 
FROM THE SOLAR FRERGY LAB AT THE UNEVERSITY OF WESCONSIN, __ 
VERSION 6.1 9/27/70 





A ree aE eT wwe wee = = @ . wee eee oe 






























SIMULATION me OsOew | 7.440€ 02 (2.590€-Q1 es ee i) | SEES 
TOLERANCES -1.000€-02 -1.000£-02 
“UNET 1 1VPE 9 WEATHER DATA READER 
PARAMETERS 6 
2-990E 0 1.CQIOE 20 2.000E 00 3.686 00 0.9 ___ 8.000€ 00 
UNIT 2 TYPE 16 RADEATION PROCESSOR, MODE 3 
ACEC ee eee - 
3.0U0E 90 36350E 02 4.050E O1 5.050t Ol 0.0 @.290E 02 2.0)0€-01 
PePuts l | 
in 2 5 : _ : _ ee ; 
0.0 
UNIT 3 TYPE 21 _ LIQUID COLLECTGR-STORAGE MUDULE _ ae oe. ee eee ee a 
PAPAMETERS 15 
e.200f 02 7.5COL 03 7.920E-O1 5.940E 03 1.Q00E 00 9.900E-O01 6.840E-01 7.000€-01 
2.120 02 Te9ODE-O1 HeOLJE OL O6260F OL  5.0290E-02 3.000E 00 Y.O)0E O1 
INPUTS) 5 
2. 1 le } Ge 1 ae 2 0, 0 
0.0 __2e6CUE OL «7 0D0E O1 2-500E 03  5.900E OL _ | oe. 
UNIT « Tyee 2y wATER TO AIR HLAT PUMP ' 
PACAMETE®SS 16 A : 
LeOVOF 99 2.500F 03 7.QJ0E O1 7. 700E O1 4.JO00E Ol 1-000E 03 8.000E 00 B.000E 00 
3.190€ O1 3.200E€ O1 L.O090E 03 oe480E 02 7.200E-01 L.100E 02 0.0 @.500€ O01 
7 INPUTS. 5 a i ee a = es Aaa. ee 
ae Me. 2 beri 5, 3 0, 0 
7.0)0F O1 2.590E 03 2.-800E Ol -2e146E 04 L-O000E 00 ; 
UNIT § TYPE 12 HEATING LUVOCENERGY/DEGREE HOUR MUDEL), MUDE 
PAFAMETERS§ A 
3.000E 00 -S-110E 02 ,7.900F O01 0.0 __ GOR nO. OSs 020) 2 00 ee 
INPUTS 4% ; . 
0. O 0, O (oh Or. O 
0.9 0.0 2-AVOE Ol 0.0 eee eee ). 
UNIT 6 TYPE 24 QUANTITY INTEGRATOR 
INPUTS 12 a ee ee | 
2. i 3, 3 Se 4% se 35 3, 7 4, 3 &e 4 &, % 
4, 6 4, 9 4,10 5, 3 | 
0.0 0.0 0.0 _ 0.0 0.0 | OvO™. _. 050 0.0 7 
a - 0.0 “6.0 °  ~ 0.0 . .}.}.}.}.»§“2e140€ 64 mak | 
. UNIT 7? TYPE 25 PRINTER(HOURLYD ao. a. oa | 
PAPAMETERS 1 ; 
1.990€ 99 
INPUTS 4 : get 
eee le scone Sy, 4e 
Rad TSTOR QuSE QLUAD | 
UNIT 8 TYPE 25 PRINTER( 26 HOUR) aa 
PARAMETERS 1 ; 
_2.4)0€ O1 a 
i a ees =. 
3s 6 6, 1 65 2 6, 3 6+ 4 6. $ 6. 6 6, 7 
_ 6+ 8 a 6102 f S.l0e eae Ogtl). = 6517 _— > 
DELU TOTRAD TUTWU WENV QTANK —~—~—\ QOUMP WORK IN QaBS 
ORES HRS QOH QAUX TOTQL 
49 nn : a ee eee 
___ EXECUTE TRNSYS ac Ie 7 06003300 ....... 
@SOMESSAGE SeLleese SUBSYSTEM TKNSYS = VERSION 86.1 05/01/78 
: ee 


TS Ne ont ba aig Te owe Ma) oa ; a at) es GD Gane “A SOS rae 8 er -_= -~< 


ey 





TEANSTENT STMULAT ION STARTING Al TIME © 0.0 
STUPPING AT TIME © 7.460€ 02 
TIMESTEP = 2.509F-OL 
: DIFFERENTIAL tWUATEIC sd ERPLK TiLEKANCE @-1 .000E-U2 




















ALGE GRATC CONVERGENCE TOLEPANLE #-1.000F-02 : 
- TIME © O10 ee a ' 
RAD TSTOR Qust QLOAD DELU TOTRAD TOTQU QENV QTANK 
0.0 Fe 000E OL 0.0 ~2.1¢0F 04 C0 a OO v3 2060 - 0.9 0.0 
; QuUMe WURK IN QABS QREJ HRS QOH QAUX TOTOL 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 ~2Z2-1%6E 04¢ 
=" TIME = 1.0000 a OSS a i ae i a 7 a a 
RAD TSTOR QUSE QLOAD 
Ope 6 T5UE O10 he 1 HOE__ 04 ——— . oe 
TIME « 2.0000 ‘ 
RAD TSTUR OUSEF —. QLUAD__ ee ee eee ht. eS ee 7 
0.0 624 66EF Ol 0.0 =~2e146E 04 
oer. TIME = SO OMA Wi We ke ee ee a 
RAD TSTUR Quse QLOAD 
{9-0 e179 Ol 0.0 =2-197E 04 
: TIME © 4.0090 . ane ae o> eon awe. * Rane . so 
‘ RAO TSTUR QUSE QLOAD 
0.0 —_—_—_ 508 D3E Cl 0.0 ee LOGE 04 ae... a a SA : 
TIME = 5.0000 
RaOd i TSTOR QusS—€ (OAD ee __ ee ee es 
0.0 5.609E OL 0.0 -2.197E 04 - 
@ 4 
- 
e 
r ? 
TIME © 742.0090 i” i. aa : : ee 
RAD TSTOR YUSE =< wl UA0 i 
0.0 Le22uel OP 0.0 whe 3 HUE 08 _ , 
TIME ® 7644.29000 | 
PAD TSTUR QUSE QUUAD i : 
0.0 Le 20kL 02 0.0 —“1L.329E 04 
TEME © 744.0300 
FAD TSTOR WUSE JL 0a DELS TOTRAD TOTQU QENV QTANA 
000 Lel7s€ 02 0.0 “Le SBDE 04 20346E 05 3.451E 04 9. 264F 06 bo 339t OS e494 06 
QDUuMP WORK IN QAIS QREYJ HRS QOH QAUX TUuTQL a 
10.0 3.399F CO -o.T6eE 06 L.O21E OF 3.908E 02 2.133€ 06 1.832E 06 -1.420€ OF 
. Gratt i. wAS CALLED 2978 TIMES | 
—— CC 2 — 2576 2° pe a ae sss Ta ii 
3 9010 ' 
= ta ~ a 6407 — AE a RE A a ye ee ee - | 
5 2378 
6 2976 
ss —— 745 
32 LETT, TI TNE eR mt 
SCOCSSCHESCSSEHTS SSCS SSS SESS OCEHEHOESESE 


¢ END WF RUN 
SCOSCHCH SHO SES ESOS SES EEE EO HEOEESEOS 





TIME © 21:32:50,802 


_ ELAPSE TIME = 186.88 SEC. 


| ee et o«.a 


ee eee 
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